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We present the control ol the high- dimensional chaos, with possibly a large number of positive 
Lyapunov-exponents, of unknown time-delay systems to an arbitrary goal dynamics. We give an 
existence-and-uniqueness theorem for the control force. In the case of an unknown system, a formula 
to compute a model-based control force is derived. We give an example by demonstrating the control 
of the Mackey- Glass system towards a fixed point and a Rossler-dynamics. 

£>y We treat the control problem for nonlinear time-delay systems. The aim of the control is to bring the 
system from a given dynamical state to an arbitrary goal dynamics. We give the conditions under which 
the control force uniquely exists. In the case of an unknown nonlinear delay system, we introduce a 
procedure to compute the control force with the help of a previously published identification procedure. 
The results are also applicable in the case of the very high-dimensional chaotic motion, with possibly a 
large number of positive Lyapunov exponents, typically observed in nonlinear delay systems. Stability 
towards noise and towards an imperfect modeling is treated in linear response. Additionally, we 
■ discuss the control of unstable periodic orbits with the help of a minimization. A numerical example, 
the control of the Mackey-Glass system from a very high-dimensional chaotic state towards a stationary 
state or a non-periodic time- varying state (the low-dimensional chaotic motion observed in the Rossler 
system), is presented in detail. 

o\ : 

I. INTRODUCTION 

C ■ 

The possibility of controlling unstable periodic orbits (UPOs) of chaotic systems with the help of suitably chosen 
external perturbances, which vanish in the case of a successful control, initiated a vast research activity. The OGY- 
method [1] and the Pyragas-method [2] attracted most attention, so far, and have been successfully applied to various 
experimental situations. While the OGY-method applies a local modeling, the Pyragas-method does not apply any 
modeling and is restricted by the value of the imaginary part of the Floquet-exponents of the unstable periodic orbit 
r* ■ to be controlled [3]. One of the aims of the control of UPOs has been to enable the switching between different 
states with the help of small perturbations only. In low-dimensional chaotic states (with only one positive Lyapunov 
exponent), though, the basic frequencies of most of the UPOs are nearly identical and a lot of UPOs are therefore 
quite similar. It is the advantage of hyperchaotic chaotic states (with more than one positive Lyapunov exponent) 
to provide for a large variety of UPOs with different basic frequencies. For this reason, the control of UPOs of 
hyperchaotic chaos is desirable. The successful control of UPOs of high-dimensional chaos with either a single control 
force [4], or a multitude of control forces [5] has been reported. More specifically, the control of of UPOs of time-delay 
systems also in hyperchaotic regimes has been achieved with a single control force only [2,6]. Despite that remarkable 
success in specific cases, it remains an open problem if all hyperchaotic states can be controlled by a single control 
force, and what actually determines the number of necessary control variables. In other words - under which necessary 
and sufficient conditions does a (scalar or multi-variate) force, which controls a hyperchaotic state towards a UPO 
exist? For the latter type of control we use the term 'control of UPOs' throughout this paper. 

On the other hand, we use the term 'control' to describe the more general task to force a system to an arbitrary, 
predefined goal dynamics, where the control force is not vanishing. The control to a predefined goal dynamics is also 
the more general, as well as the practically more relevant case, since most technical processes, such as manufacturing, 
chemical engineering, etc. rely on this type of control. It has been treated in engineering and applied mathematics 
for a long time. While for the control of linear systems a sufficient understanding has been achieved, the control of 
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nonlinear systems, in general, remains an open problem [7] . More specifically, the control of linear time-delay systems 
were summarized by Marshall [8]. Along these lines, the control of chaotic states in ordinary differential equations 
has been demonstrated [9]. In the case of an unknown time-evolution equation, though, the method has to rely on 
a modeling, i.e. with the help of Takens-like techniques, and is therefore restricted to low-dimensional chaotic states 
in typical situations. The purpose of this paper is to show that the very high-dimensional chaotic states [10] of a 
possibly unknown time-delay system with TV components can be controlled to an arbitrary goal dynamics, with A 
control forces, for the following reasons: (1) a global model of an unknown, TV-component time-delay system can 
be obtained with the help of A time series only, as shown in [11]- [13]. (2) The global model allows to compensate 
the instability-inducing effects of the time delay by suitably chosen A control forces as will be shown in this paper. 
The paper is organized as follows: At first we proove under which conditions a unique control force exists, which 
allows for a stable control of a time-delay system. This result is then used to derive a control scheme for the control 
of a time-delay system with unknown equations where a global model is achieved with an adequate identification 
procedure. Finally, we show the result of a computer experiment. 

II. THE CONTROL OF UNKNOWN DELAY SYSTEMS 

Suppose the dynamics of the system to be controlled is determined by a A-component, non-autonomous, time- 
delayed differential equation 

V = H(y,y To ,F(t)), y TQ =y(t-r Q ), (1) 

where y(t) G T>\ C 1Z N is the range of the dynamics, F(t) 6 D 2 C 7Z N is the range of the external force. Suppose the 
function H : 1Z 3N — > 1Z N is a continuous and invertible in T>2 with respect to the third argument. The control problem 
consists in proving the existence of a force F(t), which makes a predefined goal dynamics z{t) a stable solution of eq. 

(1) . In practice, global stability is desirable, since that excludes coexisting dynamical states. This can be achieved 
by imposing an additional condition on the dynamics y(t) in the form of a first-order, ordinary, non-autonomous 
differential equation 

y = g[z(t)]{y), y(0) = yo, (2) 

with g[z](y) G T>2 for y(t) G T>\. The function g : 1Z N — > 1Z N has a functional dependence on the goal dynamics 
z(t). We require z(t) to be a globally stable solution of eq. (2): i = g[z](z), with y G 2?i — ► z, for t —> oo. To 
ensure the stability of eq. (2) towards a small amount of noise (linear response regime), we also require local stability: 
5ft(Ai(z(i))) < 0,Vi, i = 1, A, where \i(z(t)) are the time-dependent eigenvalues of the time-dependent Jacobian of 
system (2). Equating the two requirements upon the dynamics we get a condition for the control force F(t), 

g[z}(y) = H(y,y T0 ,F(t)). (3) 

Since we require invertibility of the function H, the control force F(t) uniquely exists for all values of the delay time 
To and can be determined via 

F(t)=H- 1 ( y ,y ro ,g[ Z }(y)). (4) 

The control force cancels the instability-inducing effects of the time-delay. The dynamics is only determined by eq. 

(2) , and inherits its stability properties. Therefore, the goal dynamics z(t) is globally stable, and stable in linear 
response to a small amount of additional noise. But, the stability range might decrease for an increasing value of To, 
leading to an increased sensitivity towards noise for an increasing tq. 

From the practical point of view it is essential to discuss the control of delay systems, whose time-delay differential 
equations are not known, but whose dynamics y(t) is accessible to measurement and to an external driving via the 
control force F(t). For the ease of presentation, we focus on scalar time-delay systems with a single, discrete delay 
time To- We emphasize, though, that the presented ideas are also applicable in the case of multi-component time- 
delay systems, or systems with multiple delay times or a distribution of delay times. The central task is to find a 
model equation for the dynamics: F(t) = H^ 1 (y,y Tr ,y), where the function H~ x and the delay time r r are unknown 
and have to be determined from the dynamics. It has been shown, recently, that a global model of autonomous 
scalar time-delay systems can be obtained by measuring a scalar component only, even if the dynamics is very high- 
dimensional chaotic [11]- [13]. The method relies on the detection of nonlinear correlations of the observable, its time 
derivative and its time-delayed value. For the problem under consideration, a slight modification of this procedure 
is appropriate. The basic idea is to disturb the system via a time-dependent external force Fit) G Vp and measure 
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the systems reaction y(t) G T>i,y(t) € T>' x . Then, one has to find the delay time r r , for which nonlinear correlations 
in the form F(t) = H^ 1 (y,y T ,y) exist. Adequate measures to detect nonlinear correlations are ,i.e., the filling factor 
[11,12], the maximal correlation [13], or the forecast error [14]. The function H^ 1 can be determined with the help of 
adequate fitting procedures. The control force is then given by 

F(t)=H- 1 ( y ,y T0 ,g[z}(y)), (5) 

with (y,y T ) £ and g[z)(y) G T>[ if H r is invertible. There are an infinite number of functions g, which guarantee 
the existence of the control force. In practice, g can be chosen to meet the purpose best, a very simple choice will 
be: g[z](y) = —1/T(y — (z + Ti)), with a single control parameter T. If the function H r is not invertible in the 
required range, there exists no control force leading to a stable control, with the chosen g. Possibly, control can be 
achieved, then, by changing g or the action of the control upon the dynamics (via H^ 1 ). If the function H r is locally 
non-unique invertible in the required range, the control force exists in several regimes and we can define a piecewise 
continuous control force F(t) similar to eq. (5). In this case special care has to be taken for the values (y,y To ) f° r 
which the derivative of H r with respect to the third argument vanishes. 

Since, in general, the model H^ 1 slightly deviates from the actual dynamics H^ 1 = H^ 1 — eH _1 , with e being small, 
it is crucial to discuss the stability of the control towards the deviation eH -1 . For the ease of presentation, we restrict 
our discussion to the case of an additive action of the control force. Then, the dynamics of the system with control 
is determined by the non-autonomuous time-delay equation y = eh(y, y T ) + g[z](y). Calculating the deviations A(t), 
with y = z + eA in linear response, we end up with the differential equation: A = ^j(z, z)A + h(y, y T ). As expected, 
the quality of the control is directly related to the quality of the model. Deviations from the goal dynamics are 
proportional to deviations of the model as long as the linear-response approximations are valid. If not, the deviations 
of the model might lead to a failure of the control. Since we require the goal dynamics z(t) to be a globally stable 
solution of (2), and H~ x is a global model of the delayed system, the goal dynamics will, in typical cases, also be a 
globally stable solution of (1) with the control force (5). This is a huge advantage compared to a control, which only 
uses a local modeling, where the control might fail because of coexisting attractors; a situation typically encountered 
in high-dimensional chaotic systems. 



III. THE CONTROL OF THE MACKEY-GLASS SYSTEM 



In order to demonstrate the control, we apply the scheme to the Mackey- Glass system, y = f(y Ta ) — y+F(t), f(y To ) = 
tt^w , m a computer experiment for several values of the delay time r and several goal dynamics z(t). The model 



has been proposed to account for the observed large amplitude oscillations of the number of circulating white blood 
cells of patients suffering from chronic granulocyctic leukemia (first Rcf. of [19]). We assume the action of the control 
force F(t) to be additive and to be known. Therefore, it is sufficient to model the system from a time series of 
the uncontrolled system and subsequently apply the model to compute the control force F(t) via eq. (5). Because 
of the demonstrational character of the example we skip most of the details of the identification in order to keep 
the presentation short. At first, we computed a trajectory of the uncontrolled system (F = 0) for r = 80.00 with 
a Runge-Kutta algorithm of fourth order (transient time: 100 * tq; length of the trajectory: 100 * tq with 10.000 
data points). The range has been T>i = [0.20; 2.02]. The uncontrolled system exhibits a high-dimensional chaotic 
state. To estimate the delay time we used a filling factor analysis [11]. The r-dependent filling factor is shown in 
Fig. 1(a), where the delay time is indicated by a sharp local minimum. The data allow for a non-erronous estimate 
r r = 80.00 ± 0.01. In the next step, the function h r (y,y Tr ) — ,f r (y Tr ) — y has to be determined. For the fit of the 
function f r we used a rational function with a polynomial of sixth order in the nominator and a polynomial of second 
order in the denominator. The quadratic error has been 1.59 * 10~ 4 . In Fig. 1(b) the difference of the two functions 
/ = / - f r is shown. 

We choose the requirement upon the dynamics (see eq. (2)), g, in a simple linear form, 

y = g[z](y) = -l/T( y -(z + Tz)), (6) 

with a single control parameter T. The goal dynamics z(t) is a globally stable solution of (6) for y E 1Z. Therefore, 
we arrive at a non-autonomous time-delay equation for the controlled system 

y = f{y Ta )-^(y-(z + Tz)). (7) 

At first we wish to discuss the control of a fixed point, z(t) — za,z — 0, since in this case the action of the control is 
most obvious. In Fig. 2 we show the control of a chaotic state (r = 80.0) to the fixed point z = 1.5. Note that for 
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t = 80.0 the dynamical state of the uncontrolled system is very high-dimensional with 49 positive Lyapunov-exponcnts 
and a Kaplan- Yorke Dimension of 86. The deviations from the fixed point, Sy := y — z, are determined by: 



Sy = f(5y T0 + z ) - y 6 V- 



(8) 



The fixed point 5yo of (8) is given by 5yo = Tf(Sy + z ). Therefore, the magnitude of the deviations from the 
fixed point z decrease for T —> as shown in Fig. 3(a). According to eq. (8) the stability of the fixed point can 
be guaranteed for a small enough T as illustrated in Fig. 3(b). For T « 2.0, the fixed point yo undergoes a Hopf 
bifurcation to an oscillating state. 

In Fig. 4 we present the results of the successful control of the Mackey-Glass system from a low-dimensional chaotic 
state (to = 10) to a chaotic solution of the Rossler system, z(t) = xi(t), where (x\ = 0.2x\ — X2]X2 = x\ — x^;x^, = 
e + Ax%(x2 — 2)). In the case of a control towards a non-stationary goal dynamics z(t), the quality of the control is 
affected by the time scale of the goal dynamics. Loosely speaking, the parameters of the control has to be adjusted 
such that the control is faster, than the a the fastest time scale of the goal dynamics. This will be presented in more 
detail elsewhere. 



Let us shortly comment on the control of UPOs with a vanishing control force, which attracted most attention so 
far in nonlinear dynamics. The results of the paper will be also applicable in this case, if we choose the goal dynamics 
to be: z — z^ T \ where z^ T ' is a periodic orbit of the system of period T, which can be reconstructed from time series 
with well-known techniques [15]. To state our result clearly, all UPOs of a scalar time-delay system, even the ones 
which are embedded in a very high-dimensional chaotic attractor are accessible to the control with a single control 
variable only. Furthermore, we would like to emphasize that the above arguments apply to any solution of a time-delay 
system, indicating the possibility of controlling to non-stationary and non-periodic (i.e. transients, high-dimensional 
chaos) trajectories with a vanishing control force. The advantage of previously published methods for the control of 
UPOs definitely is some sort of adaptive recognition of UPOs requiring only a minimum of prior information about 
the system. In the framework of our control method, a control of UPOs of time-delay systems can be achieved by 
requiring a suitably chosen functional upon the control force F[z^](t) to be minimal under the variation of the goal 



where the functional is chosen to be an integration over the magnitude of the control force on an interval [t, t — D]. 
Technically, the minization can be performed using established nonlinear optimization procedures. 

To the authors knowledge there is no mathematical statement linking the number of necessary control variables 
to the invariant quantities of the controlled dynamics (such as attractor dimensions or number of positive Lyapunov 
exponents). It is still an open question how many control variables are necessary to ensure the unique existence of a 
control force, which gains the control of hyperchaotic chaotic states. In some cases it has been shown numerically that 
a single control force is sufficient [2,4,6]; in other cases there is evidence that more than one control force is required 
[5]. It seems to be a sufficient condition, as argued by several authors, to apply a number of control variables equal to 
the number of positive Lyapunov-exponents. Besides that in practical purposes this turns out to be unhandy, we show 
in this paper that the high-dimensional chaos of a iV-component time-delay system can always be controlled with 
a iV-component control force. Indeed, the number of positive Lyapunov-exponents is, in general, not a meaningful 
quantity in this context. We believe that the number of "localized nolinearities" as introduced in [16], which for 
time-delay systems corresponds to the number of components, might be a more appropriate concept to determine the 
number of control variables. 

Since time delays are common phenomena in nature, as well as in industrial processes, we find a wide range 
of applications for the control of time-delay systems stemming from such different fields as hydrodynamics [17], 
laser physics [18], physiology [19], and biology [20]. For instance, we envisage the possibility to suppress the large- 
amplitude oscillations of the high- Reynolds number turbulence observed in the Raylcigh-Benard system [17], with the 
time-varying temperature gradient as the control force. The same argument applies for the control of a confined jet 
[17] with the pressure difference as a control force. Equally, the control of the high-dimensional chaotic states of a 
nonlinear ring resonator [18], described by the Ikeda model, with the laser pump power as the control force, seems to 
be possible with our control method. Furthermore, in implementing the control in an adaptive manner, a control for 
non-stationary time-delay systems is envisaged. 



IV. OUTLOOK AND CONCLUSIONS 



dynamics z^ T \ i 



1. c. 




(9) 
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Finally, we summarize the conditons for the controllability of time-delay systems: (1) To control a TV-component 
time-delay system, it is sufficient to apply N control forces. (2) The function H(y,y T0 , F) has at least to be locally 
invertible with respect to F in the required range. If the underlying time-evolution of the delay system is not known, 
one has to rely on identification procedures, which additionally require: (1) The measurement of TV time series. (2) 
appropriate ranges for the identification. The control is stable against a small amout of noise, as well as structurally 
stable against small deviations of the model. 

The author acknowledges useful discussions with W. Just, H. Kantz, A. Kittel, Th. Meyer, J. Parisi, and A. Politi. 
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FIGURE CAPTIONS 



Fig. 1: (a) Filling factor under variation of the delay time r for 1000 equally sized cubes in a three-dimensional space, 
(b) The deviation of the fit f r from the function /. The dotted lines indicate the range T>\. 

Fig. 2: Control of the Mackey- Glass system for tq = 80 to the fixed point zq = 1.5. Starting from a typical initial 
condition on the attractor (transient time: 1000), the control has been switched on at t = 200 (control parameter 
T = 1) as indicated by vertical lines: (a) Time series y(t); (b) Control force F(t). 

Fig. 3: Control of the Mackey-Glass system under variation of T for r = 10 (solid circles), r = 30 (open circles), r = 50 
(stars), and r = 80 (crosses): (a) deviation of the fixed point yo from the goal dynamics z ; (b) standard- 
deviation of the controlled dynamics. 

Fig. 4: (a) Control of the Mackey-Glass system with to = 10 to a variable of the Rossler system (not shown). The 
parameter of the control has been T = 0.10. (b) Control force F(t). As indicated by lines the control has been 
switched on at t — 200. 
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